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We investigate the broadening of the bound states at an 
interface of an unconventional superconductor by bulk impu- 
rity scattering. We use the quasiclassical theory and include 
impurity scattering in the Born and in the unitarity limit. 
The broadening of bound states due to unitary scatterers is 
shown to be substantially weaker than in the Born limit. We 
study various model geometries and calculate the temperature 
dependence of the Josephson critical current in the presence 
of these impurity-broadened bound states. 



I. INTRODUCTION 

After a lot of discussion about the nature of the order 
parameter in the high-Tc materials, a number of exper- 
iments has made it virtually certain that YBa2Cu307_z 
(YBCO) is a superconductor whose order parameter has 
a d-wave symmetry. This was established by phase- 
sensitive experiments [|]J^] that investigated SQUIDs 
or ring structures containing junctions either between 
YBCO and ordinary superconductors, or between differ- 
ent domains of YBCO. 

In unconventional superconductors like those with a d- 
wave symmetry, the order parameter is sensitive to scat- 
tering from non-magnetic impurities and surface rough- 
ness. The order parameter also has a non-trivial struc- 
ture close to surfaces and interfaces. Surfaces and in- 
terfaces can be pair-breaking, i.e., the order parameter 
is suppressed on a length scale given by the coherence 
length. Andreev reflection processes of quasi-particle ex- 
citations from the spatial profile of the order parameter, 
combined with conventional (and/or Andreev) reflection 
from the surface (interface) can result under certain con- 
ditions in surface (interface) bound states. In particular, 
it has been shown || that an interface separating two 
d-wave superconductors with different signs of the order 
parameter (in a given fc-direction) will always support 
a bound state at the Fermi energy. A mathematically 
equivalent situation arises if we consider a specularly re- 
flecting wall and the sign of the order parameter is differ- 
ent for incoming and outgoing quasiparticles, i.e., there 
will also be a bound state. 

These bound states lead to a qualitative change of the 
Josephson current through a tunnel junction as has been 
recently found in It was shown that the temper- 



ature dependence of the critical current is dramatically 
different from the standard Ambegaokar-Baratoff predic- 
tion. The critical current may increase substantially at 
low temperatures or even change sign, i.e., the junc- 
tion may change its nature from an ordinary junction 
to a 7r-junction characterized by a current-phase relation 
I = Icsm((j) + 7r) = —Ic s'm(<j)). The authors of Ref. 
also considered the influence of surface roughness on 
this phenomenon and showed that it gets weaker because 
of a broadening of the bound states. In the present pa- 
per, we want to discuss how scattering from bulk impuri- 
ties changes the bound states and, as a consequence, the 
Josephson critical current. 

We use the quasiclassical formalism of superconductiv- 
ity to obtain our results. In Section || we briefly describe 
the formalism. We include bulk impurities in the stan- 
dard way by using an impurity self-energy. There are 
different models for impurity scattering, and we concen- 
trate on two limits, viz., the limits of Born scattering 
(weak scattering, scattering phase shift Sq <C 1) and uni- 
tary scattering (strong scattering, <5 — > 7r/2). 

In Section [II we calculate the angle-resolved density of 
states and the Josephson critical current in the presence 
of impurity-broadened bound states. We find by numer- 
ical calculation that bulk impurities will cut off the zero- 
temperature divergencies in the critical current predicted 
for clean systems. We also develop an analytical under- 
standing of these results and give analytical expressions 
for small scattering rates. Surprisingly, we find that the 
bound states are more sensitive to Born scatterers than 
to unitary scatterers at a given scattering rate. 



II. FORMALISM AND MODEL ASSUMPTIONS 

In this paper, we consider the system shown in Fig. [l], 
consisting of a junction between two d x 2_ y 2 superconduc- 
tors. The order parameter on side i, i = L, R is rotated 
by Ui with respect to the surface normal. The junction 
is assumed to be weakly transparent and will be mod- 
eled as a tunnel junction. We assume that the supercon- 
ductors contain a small concentration of impurities, and 
we would like to study the influence of these impurities 
on the quasiparticle bound states formed at the junction 
II- 

In unconventional superconductors, the order parame- 
ter is spatially inhomogeneous close to obstacles like sur- 
faces or interfaces. The quasiclassical formalism of su- 
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FIG. 1. Geometry of the model system: a tunnel junction 
between two misoriented d x 2_ y 2 superconductors 



perconductivity |>|-|8| is ideally suited to calculating the 
structure of the order parameter in a self-consistent way. 
Its central object, the energy-integrated Green's matrix 
function g(R, pp;e n ) fulfills the Eilenberger equation (we 
have put h = 1) 



[i£„T 3 - A(R,p F ) - E(fl, e n ),g(pF,£ n )] 

+iv F V R g(p F ,e n ) = 



(1) 



Here e n = -kT(2ti + 1) are the Matsubara frequencies, 
S(J?, s n ) is the impurity self-energy and f, are the Pauli 
matrices. 

The Green's function has to obey the normalization 
condition 



9 (PF,£n) =-TT\ 



(2) 



and A(pf) an d g(R,p F ;e n ) are related by the self- 
consistency relation (we have put ks = 1) 



A(J2,pf) - Tj2(V(PF,p' F )f(R,p' F ;e n )) 1 



(3) 



Here (...) denotes averaging over the Fermi surface, 
which we assume to be cylindrical. The model pairing 
interaction is defined by 



V(p F ,p' F ) — Acos(20— 2cti)cos(20' — 2oti 



(4) 



where i = L,R, and <fi is the azimuthal angle between p F 
and the surface normal h. 

The self-energy X describes impurity scattering, and 
we consider two models, viz., weak scattering (scatter- 
ing phase shift Sq <C 1), which will be treated in the 
Born approximation, and unitary scattering (<5o — > 7r/2). 
For an isotropic point-like impurity potential, the off- 
diagonal components of the self-energy vanish for order 
parameters transforming according to non-trivial repre- 
sentations of the point symmetry group of the crystal. 
This is the case for the d-wave order parameter studied 
here. The self-energy is then characterized by a single 
scalar function. In the Born limit the self-energy is given 
by 



X(i2; e r , 



(g(R,PF;s n )) 



1 
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(5) 



whereas in the unitary limit 
S(J2; E n ) = r t , 



(g(R,p F ;e n )) 



T u = 



ttN 



(6) 



Here, rit is the concentration of impurities, and Nq is the 
(normal) density of states at the Fermi energy. 

In the framework of the quasiclassical formalism, inter- 
faces are taken into account by Zaitsev's boundary condi- 
tions (^] . The properties of the barrier are characterized 
by the transmission probability D(p F ). 

In the limit of zero transparency (illustrated in 
Fig. |^) the boundary conditions reduce to g(pFm,0+) = 
g(pFout,0+)- We will assume specular reflection. 




FIG. 2. Geometry of incoming and outgoing quasiparticle 
trajectories for an intransparent interface. 



We are interested in the Josephson current through the 
system shown in Fig. In first order in the transparency 
this can be expressed as 

3c = - e ^Y,( v ?Min)D(pk n ) 

[f L (PF ln ^n)f + ' R (p§ out ;S n ) 

+ f + ' L U in ;e n )f R (p^ out ;e n )}) pkn , (7) 

where the Green's functions have to be evaluated at the 
interface for real order parameters, as if they would not 
have complex phases. In the results, we will eliminate 
the transparency and express it through the normal-state 
resistance of the junction which is given by (again in first 
order of the transparency) 



R N ' = e 2 A< (v^ ± {p^ in )D{p^ n )) pL 



(8) 



In the calculation, we will assume the following model 
directional dependence of the transparency: D((f>) — 
D cos(0) 2 . 



III. RESULTS FOR BORN AND UNITARY LIMIT 

A surface of an unconventional superconductor may 
support a bound state at the Fermi energy |||],[l0]] . We 
reproduce this phenomenon by a self-consistent solution 
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of a real-time version of Eqs. (0) - (^) taking into account 
impurity scattering. Figure a shows the angle-resolved 
local density of states (taken at the surface) obtained in 
the Born limit for the geometry shown in Fig. |[ The 
bound states present in a clean system can be seen to be 
broadened by impurity scattering. 
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FIG. 3. Local density of states calculated for a misorien- 
tation angle of a — 15° at a scattering rate of Tt, = 0.1Tc,o 
in Born approximation. The temperature is T = 0.01Xb,o. 
Inset: The same for a low scattering rate Tb = 0.0017b,o- 



In order to estimate the dependence of the height of 
the zero-energy peak in the density of states on the relax- 
ation time in the Born limit, we proceed analogously to 
the analytical consideration developed in jyj, based on 
the Eilenberger equations. In particular, the residue of 
the Green's function for midgap states were analytically 
calculated there in the clean limit. Due to the presence 
of impurities the pole of the retarded Green's function 
at zero energy moves into the complex plane, although 
its midgap value at the surface g(x = 0,pf;e = 0) = 
g(pp,0) is still large in the case of a sufficiently large 
relaxation time. Under this condition we obtain the fol- 
lowing integral equation for g(pp,0): 



S(ftn,F,0) X 



J dx\( g(p' F , 0) exp I - J 



dx 



,2|A(p' F ,a/ 




i xp I — / fir 





,2|A(p iriiF ,:E' 




i>xp j — J dx 



-2tt 2 \v x \t . (9) 



The order parameter can be factorized in the form 
A(pp,x) — Aoi/j(pf,x). Here ip(pp,x) is a dimension- 
less normalized function of the momentum direction and 



the distance from the surface (\i/j(pf,x)\ < 1), while A 
is the maximum value of the bulk order parameter de- 
pending upon temperature and impurity concentration. 
The spatial dependence of ip(pp,x), reduces to a depen- 
dence upon dimensionless coordinate X — A^x/vp. In- 
troducing the dimensionless variables X, X' into Eq. (||), 
we see that the propagator has the form <?(p F ,0) — 
—iy/rAoG(pF, 0), where G(pf, 0) does not contain r and 
Aq. Hence, the height of the peak in the density of states 
is proportional to \/tAq. By considering a generalization 
of the integral equation Eq. (^|) to nonzero (although suf- 
ficiently small) values of the energy, one can show ana- 
lytically that the shift of the pole position from its zero 
value is proportional to y/ Aq/t. Thus, introducing the 
dimensionless quantity fi„ = e„ a/t/Ao, we can repre- 
sent the contribution of the midgap states to the Green's 
function as g(x — 0,pf! £«) = — W T ^aG(p-p, fl n ). 

The relative strength of the influence of impurities on 
bound states in the Born and in the unitarity limits can 
be understood qualitatively by looking at Eqs. (^) and 
(||). In the absence of bound states, the Green's func- 
tion g(R,pF]£ n ) is usually quite small (compared to the 
normal-state value) for sufficiently small e n . In this case, 
according to Eqs. (||) and @, the self-energy function 
for unitary scatterers can be significantly greater than 
the one in the Born limit. By contrast, if there are 
bound states on (or quite close to) the Fermi surface, then 
the corresponding large pole-like term in the expression 
for g(R, Pf',£u) essentially rises with decreasing tempera- 
ture. This leads to the inverse situation, that is to small 
values of the self-energy for unitary scatterers as com- 
pared to the Born limit for the same values of the scatter- 
ing rates. An analogous conclusion can be drawn for the 
retarded propagator and the self-energy function taken at 
energies close to some bound state, even if it is not at the 
Fermi surface. Since the pole-like term decreases with in- 
creasing r, the above consideration, in general, does not 
work for sufficiently large values of Tb, T u . Our numerical 
calculations justify the above conclusion: in the presence 
of unitary scatterers, for sufficiently small values of T u , 
the bound states are broadened much more weakly than 
in the Born limit. In contrast to the Born limit, the 
dependence of the propagator on the parameter Ao/r„ 
does not reduce to a power-law behavior in the unitar- 
ity limit, so that simple scaling estimates are not fruitful 
in this limit. However, some rough qualitative estimates 
can be done in the unitarity limit as well, in particular, 
by considering the simplest model <i-wave order param- 
eter, whose momentum direction dependence reduces to 
A (for < 4> < 7r/2 and n < cj> < 37r/2) and to -A 
(for 7r/2 < <j> < 7r and 3n/2 < (j) < 2n). Then for small 
enough T u we find the propagator g(x = 0,pF',£n = 0) 
to be proportional to y/T u /Ao exp(AAo/r u ) with some 
numerical factor A of the order of unity. 

We will now concentrate on the temperature and impu- 
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rity dependences of the Josephson critical current. Two 
typical and experimentally relevant geometries will be 
studied: the "symmetric" junction for which qj, = an, 
and the "mirror" junction for which a/, = — cur. 

In contrast to surface roughness, bulk impurities not 
only broaden the bound states but also change the max- 
imum of the bulk pair potential and the critical temper- 
ature from their clean-system values Ao ; o and Tc,o- This 
can be seen in the lower panel of Fig. ||, where the finite 
scattering rate leads to a renormalized value of Tc- Fig- 
ure |] shows the critical current for the symmetric geome- 
try in Born approximation. The anomalous temperature 
dependence discussed in Ref. jl is still visible, but bulk 
impurity cuts off the divergence at zero temperature. 

We obtain an analytical estimate of the Josephson crit- 
ical current in the zero- temperature limit from Eq.(^|), 
taking into account the large low-energy values of the 
quantities f L ^ R \ which are associated with g the same 
way as in the clean limit fl|TT|: f(x — 0,pF',£ n ) — 
f + (x = 0,p F ;e n ) = -isignO^Aoo (pp)) g(x = 0,p F ;e n ). 
Proceeding analogously to Ref. Q for the symmetric 
(mirror) junction and introducing Q as a new integra- 
tion variable, we get 

so = ±#<v/ 2 x 



2- 



J dfi(D(j4 in )v^( P ^ n )G 2 (x = o,p£ in ,n)} p ^ . (io) 



The plus (minus) sign corresponds here to the symmetric 
(mirror) junction. 

Thus, in the Born limit and under the condition 
Aqt 3> 1 the critical current turns out to be proportional 
to A^t 1 / 2 . 

In the unitarity limit, shown in Fig. |5|, we find that the 
bound states are remarkably stable to impurity scatter- 
ing, according to our conclusion made above. For suffi- 
ciently large values of T u (when T u /Tc t o is of the order of 
unity), however, the influence can be substantial as can 
be seen in the inset to Fig. 

Figures § and @ show the corresponding results for 
the mirror junction. As in Ref. 0, the critical current 
changes sign, i.e., the junction changes character (for 
some misorientation angles) and becomes a 7r-junction 
at low temperatures. Bulk impurity scattering weakens 
this tendency. 

One simple model for a surface (or an interface) with 
roughness is a thin dirty layer containing Born impurities 
near the surface (or around the interface) Jl^-[l4|,|],|| . We 
note that the dirty layer with unitary scatterers would in- 
fluence the bound states and the Josephson critical cur- 
rent less strongly than in the Born limit (for the same 
value of the scattering rate). The problem of the most 
suitable model for surface and interface roughness for a 
given experimental situation is still open in this context. 



a, =a B 




FIG. 4. Temperature dependence of the critical current for 
the symmetric junction, i.e., at = oib. and different misori- 
entation angles a. Impurity scattering is parameterized by 
the scattering rate Ft using the Born approximation, Eq. 
Inset: Critical current for misorientation a — 45° and fixed 
temperature T = 0.005Tc,o as a function of the scattering 
rate. 
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FIG. 5. Critical current as a function of temperature for 
the symmetric junction in the unitary limit. Inset: Critical 
current for misorientation a = 45° and fixed temperature 
T = 0.01Tc,o as a function of the scattering rate. 



In conclusion, quasiparticle scattering by bulk impu- 
rities as well as surface roughness results in the broad- 
ening of surface (interface) quasiparticle bound states in 
tunnel junctions of ci-wave superconductors. In this pa- 
per, we have studied this broadening due to bulk impuri- 
ties and shown that scatterers can essentially reduce the 
height of the peak in the density of states and the low- 
temperature anomaly in the Joscphson critical current. 
We have shown that bound states are more sensitive to 
Born scatterers than to unitary scatterers at a given value 
of the scattering rate. Thus, unitary scatterers would be 
less detrimental than Born scatterers to the observabil- 
ity of the low-temperature anomalies in the Josephson 
critical current caused by bound states. 
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FIG. 6. Temperature dependence of the critical current for 
the mirror junction, i.e., ol = — or for two scattering rates 
T;, in Born approximation. 
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the mirror junction and scattering in the unitary limit. 
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